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We show that within the framework of a simple local nu- 
clear energy density functional (EDF) , one can describe accu- 
rately the one- and two-nucleon separation energies of semi- 
magic nuclei. While for the normal part of the EDF we 
use previously suggested parameterizations, for the superfiuid 
part of the EDF we use the simplest possible local form com- 
patible with known nuclear symmetries. 

PACS numbers: 21.60. Jz, 21.60.-n, 21.30.Fe 

A steady transition is taking place during the last 
several years from the mean-field description of nuclear 
properties in terms of effective forces to an energy density 
functional approach (EDF). A significant role is played in 
this transition process by the fact that an EDF approach 
has a strong theoretical underpinning [1]. The effective 
forces used to derive the EDF are nothing else but a 
vehicle, since in themselves they have no well-defined 
physical meaning. For example, the effective Skyrme 
two-particle interaction is neither a particle-hole nor a 
particle-particle interaction. The particle-hole interac- 
tion (or the Landau parameters) is defined only as the 
second order functional derivative of the total EDF with 
respect to various densities, while the particle-particle 
interaction responsible for the pairing correlations in nu- 
clei has to be supplied independently and with no logical 
connection to the Skyrme parameters. 

We shall not attempt to even mention various mean- 
field approaches suggested so far, in this respect see Refs. 
[2], but we shall concentrate instead on a single aspect 
of the nuclear EDF, namely its pairing properties. Only 
recently it became clear that a theoretically consistent 
local EDF formulation of the nuclear pairing properties 
is indeed possible [3-5]. Even though the crucial role 
of the pairing phenomena in nuclei has been established 
firmly, it is surprising to realize how poor the quality of 
our knowledge still is. Phcnomcnologically, one cannot 
unambiguously decide whether the pairing correlations in 
nuclei have a volume or/and a surface character [6-13]. 
The isospin character of the nuclear pairing correlations 
requires further clarification as well. These questions be- 
come even sharper in the language of a local EDF. 

There is also the largely practical issue of whether 
one should use a zero-range or a finite-range effective 
pairing interaction. The only reason for the introduc- 
tion of a finite-range was to resolve the formal diffi- 
culty with divergences in calculating the anomalous den- 
sities [10]. The majority of practitioners favor a much 



simpler approach, which embodies essentially the same 
physics, the introduction of an explicit energy cut-off. 
The best example are perhaps the works of the group 
[11,12], which so far is the leader in describing all known 
nuclear masses. An explicit finite-range of the pairing 
interaction tq (which can be translated into an energy 
cut-off Ec ~ h^/mrQ) and a zero-range with an ex- 
plicit energy cut-off Ec in the final analysis are equiv- 
alent. Both approaches, however, are a poor's man so- 
lution to the renormalization problem and reflect simply 
a lack of understanding of the role of high-momenta in 
the pairing channel. Neither the energy cut-off Ec nor 
the finite-range of the interaction rg carry any physical 
information and they are simply means towards getting 
rid of infinities. The argument that nuclear forces have a 
finite-range is superfluous, see Refs. [3-5], since nuclear 
pairing phenomena are manifest at small energies and 
distances of the order of the coherence length, which is 
larger than nuclear radii. 

We shall consider local nuclear EDFs only (which de- 
pend on various densities, as opposed to an explicit de- 
pendence on the full density matrix), as they proved 
overwhelmingly successful in describing normal nuclear 
properties. It is natural to expect that the same should 
apply to pairing properties. According to the general the- 
orem of Hohenberg and Kohn [1] for many fermion sys- 
tems there exists a universal EDF. Unfortunately there 
are no hints on how to derive such a functional. In the 
case of nucleons such a functional should satisfy some 
general constraints: rotational invariance, isospin invari- 
ance, time-reversal invariance and conservation of par- 
ity. Isospin symmetry is broken by Coulomb interac- 
tion, proton-neutron mass difference and charge symme- 
try breaking forces, the last two leading to rather small 
effects [13,14]. In the case of Coulomb interaction mainly 
the direct term has to be accounted for, as the exchange 
and correlation Coulomb energies seem to cancel each 
other to some extent and their combined effect together 
with the effect of charge symmetry breaking forces is rela- 
tively small and responsible mainly for such rather subtle 
effects as the Nolen-Schiffer anomaly [13,14]. We shall 
not consider here the contributions due to Coulomb ex- 
change and charge symmetry breaking energies [14] . Such 
a structure of the normal nuclear EDF i?Ar(r) would be 
complete in the absence of superfluidity. Since pairing 
correlations in nuclear systems are known to be of rather 
small amplitude, we shall consider only a superfluid EDF 
£s{f) of the following structure: 

£sir) = goir)\i^pir) + i^n{r)\^ + gi{r)\i^p{r) -Mr)\'^, 
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where fp.„(T") are the S = proton/neutron anoma- 
lous densities. There is no firm evidence of pairing in 
other partial waves except the BCS-like s-wave in ei- 
ther proton or neutron channels and the evidence for 
neutron-proton pairing is inconclusive so far. Notice that 
Ssir) is symmetric under the proton-neutron exchange. 
Wc assume that the effective couplings go.iir) might de- 
pend on position through the normal densities and that 
this dependence is consistent with expected symmetries. 
The density dependence of the effective couplings go.iC*") 
arises from two different sources. Firstly, the bare cou- 
pling constant in the pairing channel could in principle 
have some intrinsic density dependence and such depen- 
dence has been considered by various authors during the 
years [2,6-9,13]. Secondly, the renormalization of the 
pairing interaction, as described in our recent work [3-5], 
leads to position dependence as well. The equations for 
the quasi-particle wave functions Ui{r) and Vi{r) and re- 
lated quantities are 

Egs= J d^r[£N{r)+Ss{r)], 

£s{r) := -A{r)Mr) = g,ff{r)\u,{r)f , 

r [h{r) - ij]ui{r) + A{r)vi{r) = EiUi{r), 
\ A*{r)ui{r) - [h{r) - fj]vi{r) = EiVi{r), 



h{r) = -V 
1 



-V + U{r) 
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For the sake of simplicity wc do not display the spin and 
isospin variables. kp(r) is the local Fermi momentum, 
which could be either real or imaginary, while fcc(r) is 
real [3-5]. The role of the particle continuum [6,15] is 
taken into account exactly using the technique described 
in Refs. [16], the contour integration in the complex en- 
ergy plane of the Gorkov propagators for the Bogoliubov 
quasi-particles, in order to evaluate various densities. All 
calculations have been performed in coordinate represen- 
tation and all nuclei have been treated as spherical. For 
reasons we discussed in detail in Refs. [3-5], the cut-off 
energy should be chosen of the order Ec ^ 0{ep). In 
practice we found that a value Ec « 70 MeV for SLy4 
[17] and E^ « 55 MeV for FaNDF^ [13] is satisfactory 
and it ensures a convergence of the pairing field A(r) 
with a relative error of « 10~^ for density independent 
bare couplings. Note that the calculation of A(r) alone 
would require a significantly smaller Ec of order 10-15 
MeV [4,5]. The optimal value for Ec varies, depending 
on whether one uses an effective mass close to the bare 



nuclcon mass or a reduced one, as is typical with Skyrmc 
interactions. Even though this explicit cut-off energy E^ 
appears in various places, indeed, no observable shows 
any dependence on E^, when its value it is chosen appro- 
priately. Upon renormalization of the zero-range pairing 
interaction, the emerging formalism is no more compli- 
cated than a simple energy cut- off approach, with the 
only major bonus however, that there is no energy cut- 
off dependence of the results. Since the kinetic energy 
of the system is a diverging quantity of Ec and only the 
total energy is a convergent quantity [3-5,18] it is very 
important that all densities (normal and anomalous) be 
evaluated using the same energy cut-off Ec. 

We shall treat even and odd number of particles within 
the same framework and using the same EDF param- 
eterization, unlike e.g. Refs. [11,12]. The formalism for 
evaluating the Gorkov propagators for odd systems is de- 
scribed in great detail in Refs. [13,19]. For the normal 
part of the EDF we shall use either the Lyon parameteri- 
zation of the Skyrme interaction [17] and or the FaNDF*^ 
suggested by Fayans [13]. Both EDFs reproduce with 
high accuracy the infinite matter equations of state of 
Refs. [20]. 
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FIG. 1. Sn for tin and lead isotopes computed using the 
SLy4 EDF (loft) and Fayans' FaNDF" (right) with either vol- 
ume or half volume-half surface pairing. 

We shall present here results only for those nuclei for 
which we can make comparison with available recom- 
mended nuclear masses [21]. We consider at first the 
tin (38 nuclei) and lead (34 nuclei) isotope chains. We 
performed a number of calculations of these isotopes es- 
sentially from the neutron to the proton drip lines, see 
Ref. [5] for some preliminary results. For these nuclei 
we can test only the sum of the coupling constants. 
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namely g{r) = goir) + gi{r). We have considered a 
bare coupling g{r) = const, which corresponds to vol- 
ume pairing, and also g{r) = Vo[l — p{r)/pc\, with pa- 
rameters chosen to describe roiighly one half volume and 
one half surface pairing, as suggested in particular in 
Ref. [7]. One-neutron separation energies Sn and two- 
neutron separation energies 6*2,1 for tin and lead iso- 
topes were computed for constant pairing g{r) = const, 
with mean-field computed with either SLy4 interaction 
[17] (.g(r) = -250 MeV fm^) or with Fayans nor- 
mal nuclear EDF [13] {g{r) = -200MeV fm^). For 
the case of SLy4 interaction we also show results ob- 
tained for the half volume-half surface pairing model 
(Vb = -370 MeV fm^pc = 0.32 fm"^). The search 
for the appropriate values for g{r) was performed only 
among a finite set of values, e.g. in the case of volume 
pairing we considered g = —200, -225,-250,-275 and 
-300 MeV fm^ see also Ref. [5]. 
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FIG. 2. The same as in Fig. 1 but S2„. 

Table 1. Rms of 5*2^ and Sn deviations respectively 
from experiment [21] (in MeV's) for several isotope and 
isotone chains. 



Z or N 
chain 


S2n/Sn 

present 


S2n/Sn 

Ref. [11] 


S2N 

Ref. [23] 


Z = 20 


0.82/0.76 


1.02/0.92 


0.96 


Z = 28 


0.67/0.50 


0.66/0.55 


1.30 


Z = 40 


0.93/0.63 


().()6/0.63 


2.21 


Z = 50 


0.29/0.21 


0.43/0.35 


0.95 


Z = 82 


0.23/0.37 


0.58/0.53 


0.74 


N = 50 


0.37/0.26 


0.41/0.23 


NA 


N = 82 


0.43/0.31 


0.50/0.56 


NA 


N = 126 


0.42/0.23 


0.88/0.52 


NA 



The agreement between experiment and theory is par- 



FIG. 3. Sp for isotone chains N = 50, 82 and 126 and S„ 
for calcium isotopes. For S2p and S2n see online report [22]. 



ticularly good for the case of FaNDFO. We relate this 
result with the fact that the effective mass in FaNDF° is 
the bare nucleon mass, unlike the case of SLy4 EDF and 
in agreement with the global mass fit of Ref. [11,12]. It 
is notable that the agreement with experiment is equally 
good for both tin and lead isotopes with the same value 
of g, unlike Refs. [13]. Even though we used the same 
normal nuclear EDF as in Refs. [13], our agreement with 
experiment is notable superior, see Ref. [22], even though 
we parameterize the pairing interaction with one param- 
eter only vs. up to five parameters used in these papers. 

Since FaNDF*^ in conjuction with the bare pairing cou- 
pling constant g{r) = —200 MeV fm^ apparently pro- 
vides the best description in case of tin and lead isotopes, 
further calculations were performed only with this choice 
of parameters. In Fig. 3 we display the one-proton sep- 
aration energies Sp for three isotone chains (23 nuclei 
with A'' = 50, 25 nuclei with N = 82 and 14 nuclei with 
N = 126) and Sn for calcium isotopes (24 nuclei), and 
for nickel and zirconium isotopes see Ref. [22], 212 nu- 
clei in total. Since the neutron numbers for these isotone 
chains are also magic, again, we can test only the same 
combination of coupling constants g{r) = 3o('") + 5i (''')• 
As we have conjectured at the beginning of this study, 
one can indeed describe with a single value g separately 
for proton and neutron pairing correlations in both even 
and odd systems, as opposed to the treatment of Refs. 
[11,12], which slightly violates isospin invariance. In es- 
sentially all cases in which we have been able to perform 
a comparison between our results and those available in 
literature, our results were either qualitatively superior 
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or, in a few separate cases, as good as any other results. 
In Table 1 we present rms deviation from experimental 
(recommended) values [21] for the two and one nucleon 
separation energies for several isotope and isotonc chains. 
The size of each set of nuclei in a chain was given by the 
number of nuclei in Ref. [11], for which there are experi- 
mental values in the unpublished Audi and Wapstra 2001 
compilation. 

There are a number of theoretical arguments, suggest- 
ing that the pairing coupling should be density /position 
dependent, due to the coupling to surface/particle-hole 
modes, e.g. Ref. [24]. A similar line of reasoning was 
presented in the case of dilute systems [25,26] and neu- 
tron matter [27,28] for quite some time. Our results, see 
Figs. 1-2, show that Sn and S2n for tin and lead isotopes 
are not particularly sensitive to such effects. To some ex- 
tent this is not a surprise, since pairing correlations are 
"built" at distances of the order of the coherence length 
S, oc h^kp/mA [29], see Ref. [30] for a related instructive 
example. This apparent low sensitivity of Sn and <S'2jv 
to a possible density dependence of the pairing couplings, 
could in principle be profitably used to describe other ob- 
servables. From the results of Refs. [26] one might infer 
that pairing coupling constants could have a noticeable 
variation with the isospin composition of a given system, 
since the magnitude of the induced interactions changes 
dramatically as the number of fermion species varies. 
Neither our results, nor previous work has necessitated 
the introduction of such a dependence however. In our 
phenomenological approach, based on general symmetry 
arguments alone and the fact that the pairing correlations 
are relatively weak in nuclear systems, we restricted the 
form of the EDF superfluid contribution to the simplest 
one compatible with known symmetries. We were able 
to infer that pairing properties of either kind of nucleons 
can be accounted for with a single constant g = go+gi- It 
remains to be seen whether the other (non-perturbative) 
combination g' = go — gi (never considered by other au- 
thors) could ever become relevant. 
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